Pairing resonance as a normal-state spin probe in ultra-thin Al films 
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We present a quantitative analysis of the low-temperature, high parallel field pairing resonance 
in ultra-thin superconducting Al films with dimensionless conductance g S> 1. In this regime we 
derive an analytical expression for the tunneling density-of-states spectrum from which a variety 
of normal-state spin parameters can be extracted. We show that by fitting tunneling data at 
several supercritical parallel magnetic fields we can determine all of the relevant parameters that 
have traditionally been obtained via fits to tunneling data in the superconducting phase. These 
include the spin-orbit scattering rate, the antisymmetric Landau parameter C", and the orbital 
pair-breaking parameter. 

PACS numbers; 74.50.-|-r, 74.40.-l-k, 74.78.Db, 73.50.Fq 



I. INTRODUCTION 

Determining the microscopic spin parameters of 
paramagnetic metals has historically been a process 
fraught with comphcations and inaccuracies.^'- In gen- 
eral, the spin response of an interacting fermionic sys- 
tem can be modified by spin-orbit scattering processes, 
electron-phonon interactions, and/or electron-electron 
interactions.^'"* These contributions to the spin sus- 
ceptibility themselves can be affected by disorder;^ 
dimensionality,^'^ and the presence of interfaces!^ The 
two primary spin parameters for a paramagnetic system 
are the spin-orbit scattering rate and the antisymmet- 
ric I = Landau parameter G^. The latter accounts 
for the renormalization of the bare Pauli spin suscepti- 
bility due to electron-phonon and electron-electron in- 
teractions. Depending upon the sign of this parame- 
ter the effective spin moment can be larger or smaller 
than the bare electron value. In practice, the spin- 
orbit scattering rate can be obtained from the coherent 
backscattering contributions to the magneto-resistance 
of moderately disordered non-superconducting films or 
by parallel magnetic field studies of thin superconduct- 
ing films. The Fermi-liquid parameter G'^, however, is 
more difficult to determine accurately. In principle, it can 
be extracted from low-temperature measurements of the 
spin susceptibility x and the heat capacity 7 from which 
the respective corresponding density of states N{x) and 
7V(7) are obtained. The ratio of these densities of states 
is a direct measure of the many body renormalization, 
G'^ — N{'f)/N{x) — 1^2 Unfortunately, orbital contribu- 
tions to the susceptibility make it very difficult to deter- 
mine its spin component precisely in high-conductivity 
systems and phonon contributions to the specific heat 
can introduce significant systematic errors in the mea- 
surement of iV(7). In this report we address the deter- 
mination of G*^ and the spin-orbit scattering rate via the 
Pauli-limited, normal-state pairing resonance FiSiiiiiSiiH 

If a paramagnetic system has a superconducting phase 



and can be made into a thin-film form, then it is pos- 
sible to access the spin parameters through tunneling 
density-of-states (DOS) measurements. A Zeeman split- 
ting of the BCS coherence peaks can be induced by ap- 
plying a parallel magnetic field to a film of thickness 
t <ti where ^ is the superconducting coherence length. 
Tedrow and Meservey pioneered the use of supercon- 
ducting spin-resolved tunneling to directly measure both 
spin-orbit scattering rate and the Landau parameter G° 
in thin Al and Ga films near the parallel critical-field 
transition,- '^^'^^ This technique, however, cannot access 
G° well into the superconducting phase since those elec- 
trons responsible for the exchange effects are consumed 
by the formation of the condensate To circumvent this 
limitation, one needs to measure the Zeeman splittings in 
magnetic fields just below parallel critical field. However, 
one cannot completely reach the normal-state quasiparti- 
cle density in a thin film while remaining in the supercon- 
ducting phase, since the spin-paramagnetically limited 
parallel critical-field transition is first-order. Because of 
this, one must extrapolate the normal-state value of G'^ 
from data taken in the superconducting phase. Alterna- 
tively, the films can be made marginally thicker, which 
will suppress the first-order transitioufi^ or the measure- 
ments can be made at higher temperatures. But these 
strategies limit one to a very narrow range of film thick- 
nesses. Furthermore, in both cases one is constrained to 
a very narrow range of applied fields. 



Here we present a detailed analysis of the normal-state 
pairing resonance (PR) from which the spin-orbit scat- 
tering rate, orbital depairing parameter, and the Lan- 
dau parameter G" can be accurately obtained. We show 
that the technique can be used over a wide range of 
film thicknesses and resistances. Moreover, the measure- 
ments can be made in fields well above the parallel critical 
field and in fields substantially tilted away from parallel 
orientatioui^ii^ 
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II. PAIRING RESONANCE IN PARALLEL 
FIELD 

The PR is characterized, as any other resonance, by 
two quantities: its position and its width. The former is 
given byii 



E+^\ {Ez + n) , 



where 



1 + GO 



(1) 



(2) 



is the Zeeman energy renormahzed by the Fermi-Hquid 
parameter G^, hb is the Bohr magneton, and 



El-^l 



(3) 



is the Cooper-pair energy with Aq the zero-field, zero- 
temperature gap of the corresponding superconducting 
phase. 

The width of the PR depends on the effective dimen- 
sionahty of the sample and on the strength F of pair- 
breaking mechanisms other than the Zeeman splitting. If 
these are absent, a non-perturbative approach is neces- 
sary (see Ref. [ill), and for quasi- two dimensional systems 
the width is 



W2 = 



A2 



(4) 



where g = A'kKvqD is the dimcnsionless conductance with 
D the diffusion constant and the bare DOS. If <C F 
then a perturbative calculation is sufficient to accurately 
estimate the width, provided one properly takes into ac- 
count the role of the exclusion principle.^'* For instance, 
in the case of a tilted magnetic field, F is proportional 
to the perpendicular component of the field and the ex- 
clusion principle both shifts and reshapes the PR. If we 
consider the effects of spin-orbit scattering and the finite- 
thickness orbital contributions of the parallel field, then 



F 

2A;; 



= 6 



^^BH 
Ao 



(5) 



where, according to the notation commonly used to char- 
acterize the DOS in the superconducting sisXe^ 



b = 



(6) 



is proportional to the spin-orbit scattering rate 1 / Tso and 

DeH^Ao 



(7) 



is the orbital de-pairing parameter, where t is the film's 
thickness, e is the electron charge, and I is the mean-free 
path. This latter parameter quantifies the strength of the 



orbital effect of the field-- in relation to the Zeeman ef- 
fect. The Zeeman splitting is the dominant pair-breaking 
mechanism for c < 1. 

Following the procedure outlined in Ref. 18, we obtain 
the zero-temperature correction to the (spin-down) DOS 
due to the PR 



= -A(e;i?z,r) 



F2 



(8) 



where e is the energy measured from the Fermi level; the 
other quantities entering this formula have been defined 
above, see Eqs. ((J)-®. The correction for the other spin 
component is found by replacing e ^ — e in the right- 
hand side of Eq. The function 

A{e;Ez,r) = -\ ai-ctanHEz - e)/F] arctan[f]/F] 

-Karctan[(e - f7)/F] 4- arctan[(2e - Ez)/T]j 

(9) 

accounts for the exclusion principle and takes on values 
between and 2. It alters the Lorentzian shape of the 
PR, especially at energies close to the Fermi energy (i.e., 
e < E+) and, in fact, A{e = 0) = 0. We note that 
Eqs. (E])-© imply that Sv/vq < 2W2/F, which is consis- 
tent with the assumed perturbative criterion F 3> W2. 

In this work we show that Eq. ([5]) gives a quantitative 
description of the PR and that it enables us to extract 
from normal-state measurements the physical quantities 
G", b, and c. While they can be obtained from DOS mea- 
surements in the superconducting state,- this requires 
to solve a set of self-consistent equations for the order pa- 
rameter and "molecular" magnetic field together with the 
Usadel equations for the normal and anomalous Green's 
functions - a much more complicated task in compari- 
son to the simple fitting of the data that we describe in 
Sec.lTVl 



III. EXPERIMENTAL PROCEDURE 

Aluminum films were grown by e-beam deposition of 
99.999% Al stock onto fire-polished glass microscope 
slides held at 84 K. The depositions were made at a rate 
of ~ 0.1 nm/s in a typical vacuum P < 3 X 10""^ Torr. 
A series of films with thicknesses ranging from 2 to 
2.9 nm had a dimcnsionless normal-state conductance 
that ranged from g — 5.6 to 230 at 100 mK. After deposi- 
tion, the films were exposed to the atmosphere for 10-30 
min in order to allow a thin native oxide layer to form. 
Then a 9-nm thick Al counterelectrode was deposited 
onto the film with the oxide serving as the tunneling bar- 
rier. The counterelectrode had a parallel critical field of 
~2.7 T due to its relatively large thickness, which is to 
be compared with Hi,\\ ^ 6 T for the films. The junc- 
tion area was about 1 mmxl mm, while the junction re- 
sistance ranged from 10-100 kf2, depending on exposure 
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FIG. 1; Tunneling conductance at 70 mK for three super- 
critical parallel magnetic fields (solid lines). The dashed line 
is the best fit to the zero bias anomaly due to Coulomb in- 
teraction. The arrows point to the boundaries of the low- 
{\V\ < 0.2 mV) and high-bias {\V\ > 1.4 mV) regions used 
for the fitting. 



FIG. 2: Pairing resonance at 7 T (solid line) with the ZBA 
subtracted off. The dot-dashed curve is the best fit to the data 
using Eq. (|8]). The dashed curve is the best fit with a sum 
of Eq. (|8]) and a Gaussian - see the text for more details on 
the fitting procedure. The two terms of the sum are plotted 
separately as dotted curves. 



time and other factors. Only junctions with resistances 
much greater than that of the films were used. Mea- 
surements of resistance and tunneling were carried out 
on an Oxford dilution refrigerator using a standard ac 
four-probe technique. Magnetic fields of up to 9 T were 
applied using a superconducting solenoid. A mechanical 
rotator was employed to orient the sample in situ with a 
precision of ^ 0.1°. 



IV. RESULTS AND DISCUSSION 

We show in Fig. [T] the tunneling conductance mea- 
sured at 70 mK and three supercritical parallel magnetic 
fields. This particular film of dimensionless conductance 
g ~ 57 was 2.6 nm thick and had a zero-field supercon- 
ducting transition temperature Tc = 2.74 K. Common 
to the three data sets is the Coulomb zero-bias anomaly 
(ZBA),^^ which produces a logarithmic depletion in the 
DOS at high biases; the logarithm is cut off at low bias 
by temperature. In order to isolate the paramagnetic res- 
onance, we need to remove the contribution of the ZBA. 
To interpolate between the low- and high-bias parts of 
the curves (as delimited by the arrows in Fig.[l]), we find 
the best-fit curve, restricted to these regions, given by the 
sum of a background constant tunneling conductance and 
Re ^'(1/2 -I- iaV), where ^ is the digamma function and 
a a fitting parameter. The result is the dashed curve 
in Fig. [U which is then subtracted from the measured 
tunneling conductances. 

In Fig. [2] we plot with a solid line the PR at 7 T ob- 
tained as described above. As discussed in Sec. HIl its 
position and width are respectively determined by the 



Zeeman energy Ez and the pair-braking rate T, while the 
conductance g only affects the overall magnitude. Using 
Eq. ([8|), the best fit to the data is given by the dot- 
dashed curve; while the main peak is well reproduced, 
a shoulder feature at higher bias is underestimated. To 
our knowledge, there are two possible causes for this dis- 
crepancy, namely a finite bias, triplet channel anomaly,— 
similar to the Coulomb ZBA but much weaker, and finite- 
temperature effectsj^ To take into account these possible 
corrections, we add to Eq. ([8]) a Gaussian contribution; to 
reduce the number of free parameters, we require it to be 
centered at the Zeeman energy, which is where a triplet 
channel correction would be located, while the amplitude 
and width are used as fitting parameters. The best fit 
thus found is the dashed line in Fig. [21 the peaked PR 
and broad Gaussian contributions are plotted separately 
with dotted lines. 

We present in Fig. [3] two more PRs with the best-fit 
curves. The asymmetric shape of the resonance and its 
suppression near the Fermi energy are evident in the low- 
est field data. We note that fitting these data with Eq. ([5]) 
only would require us to decrease the conductance with 
increasing field, whereas we can use the same value of the 
conductance at all fields when the Gaussian correction is 
included. Moreover, the value of the Zeeman energy is 
only weakly affected by the inclusion of this correction, 
with the change in Ez smaller than our estimated rel- 
ative error of about 1%. While these two observations 
support the validity of our approach, the magnitude of 
the width parameter T turns out to be more sensitive to 
the Gaussian correction. However, its field dependence 
(see Fig. [5]) is robust, and the quantitative estimates dis- 
cussed below are in line with expectations. 
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FIG. 5: Normalized pair-breaking parameter F/Aq vs. the 
square of the reduced field. Using the linear relationship in 
Eq. (O we obtain from the best-fit line the spin-orbit scatter- 
ing rate b ~ 0.06 and the orbital effect parameter c ~ 0.02. 
As in Fig. 3] we show in the inset the data pertaining to the 
2.9 nm-thick film. 



FIG. 3; Pairing resonances measured at 8 T (top) and 6 T 
(bottom solid curve). The bottom curve is shifted down by 
0.007 for clarity. The dashed lines are best fits to the data ob- 
tained as described in the text. The asymmetry of the PR and 
its suppression near the Fermi energy are easily recognized in 
the data taken at 6 T. 




FIG. 4: Normalized Zeeman energy Ez /Ao vs. magnetic 
field H. The solid line is the best fit to Eq. ©; the slope 
is proportional to (1 -I- G°)~^, and we estimate the value of 
the Fermi-liquid parameter G° ~ 0.26. For comparison, the 
dashed line represents the expected linear relationship in the 
absence of Fermi-liquid renormalization. Inset: same plot as 
the main figure, but for a thicker film with G" ~ 0.24 (see 
text for details). 



Having detailed our fitting procedure, we now consider 
the physical quantities that can be extracted from the 
data. In Fig. [4] we plot the normalized Zeeman energy as 
a function of the applied field. By fitting the data with 
Eq. (HI) we find G" ~ 0.26; a similar estimate, G° ~ 0.24, 
is obtained for a thicker film with t = 2.9 nm, g — 230, 
and zero-field, zero-temperature gap Aq — 0.41 meV - 
see the inset of Fig. 21 We note that a better fit to 
the data in Fig. [4] could be obtained by allowing for a 
finite negative intercept; however, the large estimated 
error on the intercept makes the best-fit line compati- 
ble with the expectation that it passes through the ori- 
gin [see Eq. ([2])]. This finite intercept could be due to 
small higher-order contributions, since at the lowest field 
the parameter 2W2/T ~ 0.07 is only marginally smaller 
than 1. In support to this interpretation, we find no ev- 
idence of finite intercept for the thicker film for which 
2W2/T < 0.016. Alternatively, the intercept could be an 
additional indication, together with the shoulder feature 
mentioned above, of finite-temperature effects. We will 
further investigate this latter issue in a separate work. 

The width parameter F is plotted in Fig.[5]as a function 
of (fiBH/Ao)^ together with the best-fit Hue. According 
to Eq. ([5]), the intercept and the slope are determined 
by the spin-orbit parameter b and orbital parameter c, 
respectively. We estimate their values as 6 ~ 0.06, in 
agreement with the results in the literature, and c ~ 0.02, 
which favorably compareg^^ with the value c ~ 0.04 ex- 
trapolated from superconducting state measurements in 
marginally thick (i.e., c ~ 1) films. Repeating the analy- 
sis for the thicker film - see the inset of Fig. [5] - we find 
b ~ 0.06 and c ~ 0.04. As a further check on the validity 
of the present approach, for this film we show in Fig.[5]the 
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FIG. 6: Tunneling DOS in the normal (left, H = 5.6 T) and 
superconducting (right, H = 4 T) states at T = 70 mK for a 
2.9 nm thick film. Solid lines are experimental data; dashed 
lines have been calculated with the parameters given in the 
text. 

measured and calculated DOS in the normal and super- 
conducting states for fields of 5.6 and 4 T, respectively: 
all the main features of the superconducting DOS are 
captured by the theoretical curve^^ obtained by solving 
the Usadel and self-consistent equations"'^^ with the pa- 
rameters found via the normal-state measurements. 



In summary, we have presented a quantitative study of 
the paramagnetic pairing resonance in parallel field. We 
have derived an expression, Eq. ([5]), for the density of 
states which takes into account spin-orbit scattering, or- 
bital effect of the magnetic field, and the Pauli exclusion 
principle. The latter is responsible for the suppression 
of the resonance near the Fermi energy, see Fig. [3] and 
the left panel of Fig. [51 By fitting the PRs measured at 
different fields we have obtained the values of the Fermi- 
liquid parameter , the spin-orbit scattering rate &, and 
the orbital parameter c, thus showing that normal-state 
experiments can provide the same information usually 
extracted from the DOS of the superconducting phase. 
Since the PR affects the spin-resolved DOS at opposite 
biases, it can, in fact, be used to probe the electron-spin 
polarization in magnetic films. The present work pro- 
vides the foundation for the analysis of tunneling studies 
of itinerant magnetic systems via the PRi^ 
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